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ABSTRACT 

Kepler will monitor a sufficient number of stars that it is likely to detect 

single transits of planets with periods longer than the mission lifetime. We show 

that by combining the exquisite Kepler photometry of such transits with precise 

radial velocity observations taken over a reasonable timescale (~ 6 months) after 

the transits, and assuming circular orbits, it is possible to estimate the periods 

of these transiting planets to better than 20%, for planets with radii greater 

than that of Neptune, and the masses to within a factor of 2, for planets with 

masses larger than or about equal to the mass of Jupiter. Using a Fisher matrix 

analysis, we derive analytic estimates for the uncertainties in the velocity of the 

^ \ planet and the acceleration of the star at the time of transit, which we then use 

to derive the uncertainties for the planet mass, radius, period, semimajor axis, 

and orbital inclination. Finally, we explore the impact of orbital eccentricity on 

O ' the estimates of these quantities. 

00 ' 
O 



X 



Subject headings: methods: analytical, planetary systems, planets and satellites: 
general 



1. Introduction 

Planets that transit their stars offer us the opportunity to study the physics of planetary 
atmospheres and interiors, which may help constrain theories of planet formation. From the 
photometric light curve, we can measure the planetary radius and also the orbital inclina- 
tion, which when combined with radial velocity (RV) observations, allows us to measure the 
mass and density of the planet. Infrared observation s of planets during secondary eclipse can 



be us ed to measure the planet's thermal spectrum (ICharbonneau et al.ll2005l ; iDeming et al. 



20051 ) . and spectroscopic observations around the time of transit can constrain the com- 



position of the planet's atmosphere through transmission spectroscopy (ICharbonneau et al. 
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2002 ). Additional ly, optical observations of the secondary eclipse can probe a planet's albedo 



( IRowe et al.l 120071 ). All of the known transiting planets are Hot Jupiters or Hot Neptunes, 
which orbit so close to their parent stars that the stellar flux plays a major role in heating 
these planets. In contrast, the detection of transiting planets with longer periods (P > lyr) 
and consequently lower equilibrium temperatures, would allow us to probe a completely dif- 
ferent regime of stellar insolation, one more like that of Jupiter and Saturn whose energy 
budgets are dominated by their internal heat. 

Currently, it is not feasible to find long-period, transiting planets from the ground. 
While RV surveys have found ~ 110 planets with periods > lyr, it is unlikely that any 
of these transit their parent star, given that for a solar- type host, the transit probability 
is very small, p tI ~ 0.5%(P/yr)~ 2//3 . Thus the sample of long-period planets detected by 
RV will need to at least double before a transiting planet is expected. Given the long 
lead time necessary to detect long-period planets with RV observations, this is unlikely to 
happen in the near future. Furthermore, the expected transit times for long-period systems 
are relatively uncertain, which makes the coordination and execution of photometric follow- 
up observations difficult. Putting all this together, radial velocity surveys are clearly an 
inefficient means to search for long-period, transiting planets. 

Ground-based photometric transit surveys are similarly problematic. A long-period 
system must be monitored for a long time because typically, at least two transits must be 
observed in order to measure the period of the system. Furthermore, if there are a couple of 
days of bad weather at the wrong time, the transit event will be missed entirely, rendering 
years of data worthless. Thus transit surveys for long-period planets require many years of 
nearly uninterrupted observations. Additionally, many thousands of suitable stars must be 
monitored to find a single favorably-inclined system. 

Because of its long mission lifetime (L = 3.5 yrs), near ly continuous observations, and 
large number of target stars (N ~ 10 5 ), the Kepler satellite (IBorucki et al.ll2004j ; iBasri et al. 



20051 ) has a unique opportunity to discover long-period transiting systems. Not only will 
Kepler observe multiple transits of planets with periods up to the mission lifetime, but it is 
also likely to observe single transits of planets with periods longer than the mission lifetime. 
As the period of a system increases beyond L/2, the probability of observing more than 
one transit decreases until, for periods longer L, only one transit will ever be observed. For 
periods longer than L, the probability of seeing a single transit diminishes as P~ 5 / 3 . Even 
with only a single transit observation from Kepler, these long-period planets are invaluable. 

We show that planets with periods longer than the mission lifetime will likely be de- 
tected and can be characterized using the Kepler photometry and precise radial velocity 
observations. Furthermore, this technique can be applied to planets that will transit more 
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than once during the Kepler mission, so that targeted increases in the time-sampling rate 
can be made at the times of subsequent transits. 

In §2] we calculate the number of one- and two-transit systems Kepler is expected to 
observe. We give a general overview of how the Kepler photometry and radial velocity 
follow-up observations can be combined to characterize a planet in £j3l We discuss in £jH the 
expected uncertainties in the light curve observables, and we relate these to uncertainties in 
physical quantities, such as the period, that can be derived from the transit light curve. §5] 
describes the expected uncertainties associated with the RV curve and how they influence 
the uncertainty in the mass of the planet. $6] discusses the potential impact of eccentricity on 
the ability to characterize long-period planets. We summarize our conclusions in ^3 Details 
of the derivations are reserved for an appendix (§HD- 



The number of single-transit events Kepler can expect to find is the integral of the 
probability that a given system will be favorably inclined to produce a transit, multiplied by 
the probability that a transit will occur during the mission, convolved with a distribution 
of semimajor axes, and normalized by the expected frequency of planets and the number of 
stars that Kepler will observe, 



Here N tot is the total number of transiting systems expected, is the number of stars being 
monitored, jp tr is the probability that the system is favorably inclined to produce a transit 
(the transit probability), pi is the probability that a transit will occur during the mission, 
and f(a) is some distribution of the planet semimajor axis a, normalized to the expected 
frequency of planets. 

Assuming a circular orbit, the transit probability is simply, 



where .R* and M* are the radius and mass of the parent star, and P is the period of the 
system. 



2. Expected Number of Planets 




(1) 




(2) 



We consider both the probability of observing exactly one transit, pL,i, and the proba- 
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bility of observing exactly two transits, pL,2- These probabilities are 
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Thus there will be a range in periods from L/2 to L where it is possible to get either one 
or two transits during the mission. If we assume a mission lifetim^] of L = 3.5 years and a 
solar-type star (R* = R & , M* = M ), we can find the total probability (ptrpi) of observing 
a planet that transits and exhibits exactly one or exactly two transits as a function of period 
or semimajor axis. These probability distributions as functions of a are shown in the middle 
panel of Fig. [IJ 



For P > L, the total probability, p to t, of observing a single transit is 



Ptot = ptvpL = 0.002 
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(5) 



Thus the probability that the planet will transit and that a single transit will occur during 
the mission is generally small for periods at the mission lifetime, and beyond this is a strongly 
decreasing function of the period, oc P~ 5 / 3 or oc a~ 5 ^ 2 . 

The observed distribution of semimajor axes for planets with minimum masses m p sin % > 
0.3 Mj up is shown in the top panel of Figured], where m p is the mass of the planet, i is the 
orbital inclination, and the data are taken from the "Catalog of Nearby Exoplanets" of 
Butler et al. Note that we have included multiple-planet systems. We normalized 

the distrib ution by assuming tha t 8.5% of all stars have planets with m p sini > 0.3 Mj up and 
a < 3 AU (jCumming et al.ll2008l ). We can extrapolate this distribution to larger semimajor 
axes by observing that it is approximat ely constant as a function of log a and adopting 
the value found by lCumming et al.l (120081 ). dN/dloga = 4.3%. Using this extrapolation, we 
predict 13.4% and 16.4% of star s have a m p sin % > 0.3 Mj, ip planet within 10 AU and 20 AU, 
respectively, in agreement with ICumming et al.l (120081 ) . 



We find the number of planets expected for Kepler by convolving this normalized distri- 
bution with the probability distribution and multiplying by N* = 100, 000. We convolve the 



1 Here and throughout this paper, we use the characteristics of the Kepler mission provided by the website, 



http://kepler.nasa.gov/| 



2 The catalog is available from http:/ /exoplanets. org/planets.shtml It was accessed 2008 May 9. 
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distribution of semimajor axes with the total transit probability distribution by weighting 
the probability of observing a transit for each member of the bin to determine an overall 
probability for the bin. The bottom panel of Fig. [T] shows the resulting distribution of one- 
and two-transit systems. Integrating this distribution over all semimajor axes out to the 
limit of current RV observations predicts that Kepler will observe a total of 4.0 single-transit 
systems and 5.6 two-transit systems during the mission lifetime of L = 3.5 yrs. If we include 
the extrapolation of the observed distribution to larger semimajor axes, the integrated num- 
ber of single-transit systems increases to 5.7; the effect is modest because the probability of 
observing a transit declines rapidly with increasing semimajor axis (p 

tot CX CL J. 

Thus, Kepler is likely to detect at least handful of single-transit events. It is important 
to bear in mind that this is a lower limit. First, our estimate does not include planets with 
m p smi < 0.3 Mj up , primarily because RV surveys are substantially incomplete for low- 
mass, long-period planets. As we will show, Kepler will be able to constrain the periods of 
planets with radii as small as that of Neptune with the detection of a single t ransit. Indeed, 



microlensing surveys in dicate tha t cool, Neptu ne-m ass planets a re common (IBeaulieu et al. 



2006l : lGould et al.lbood ). Second, liarnesl (ho07h andEurfl (hoosh showed that a distribution 



in eccentricity increases the probability that a planet will transit its parent star. Finally, if 
the Kepler survey lifetime is extended, the number of detections of long-period transiting 
planets will also increase. 

One might assume that a single-transit event must simply be discarded from the sample 
because the event is not confirmed by a second transit, and the period cannot be constrained 
by the time between successive transits. Although we do not expect to see many of them, 
long-period transiting planets are precious, as the information we could potentially gain by 
observing these systems when they transit in the future would allow us to greatly enhance 
our understanding of the physical properties of outer giant planets, and in particular would 
allow us to compare them directly to our own solar system giants. Thus, these single transit 
events should be saved, if at all possible. 



3. Characterizing a Planet 

With a few basic assumptions, Kepler photometry of a transiting planet can provide a 
period for the system that can be combined with precise radial velocity measurements to 
estimate a mass for the planet. Here we briefly sketch the basics of how these properties 
can be estimated and what their uncertainties are, and we provide details in the following 
sections. We assume the planet is much less massive than the star, m p <C M*, and circular 
orbits. We consider the impact of eccentricity on our conclusions in §HJ 
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Fig. 1. — {top) The observed distribution of semimajor axes for planets with minimum masses 
m p smi > 0.3 M Jup . The shaded region indicates a > 10 a5 AU ~ 3AU, where we expect 
the sample to suffer from incompleteness. The dashed line indicates our extrapolation of 



the di stribution, found by assuming the constant distribution in log a from ICumming et al. 



(120081 ) . (middle) The total transit probability. The solid line shows the probability that a 
planet transits and exhibits a single transit during the Kepler mission lifetime. The dotted 
line shows the same probability for two transits, (bottom) The expected number of one 
and two-transit events. The probability distribution has been convolved with the observed 
fraction of planets and multiplied by the number of stars Kepler will observe (100,000) to 
give the number of single transit systems (solid line) and two transit systems (dotted line) 
that can be expected. The dashed line shows the convolution of the single transit probability 
distribution with the extrapolation to larger semimajor axes. 
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Assuming no limb-darkening and assuming the out-of-transit flux is known perfecthjfl, a 
transit can be characterized by four observables, namely the fractional depth of the transit 
S, the full-width half-maximum of the transit T, the ingress/egress duration r, and the time 
of the center of transit t c . These can be combined to estimate the instantaneous velocity of 
the planet at the time of transit, 

v tI , p = 2rJ^) . (6) 



For a circular orbit, v tT , P can be used to estimate the period of the planet 

„ 2na Gn 2 /Tr\ 3/2 

■P*\-r?) > ( 7 ) 



v tr,p 3 \ ^/S J 

where we have employed Kepler's third law assuming that the planet's mass is much smaller 
than the star's mass. Thus, if the stellar density p+ is known, then the planet's p eriod can be 



estimated from photometry of a single transit (ISeager fc Mallen-Ornelad 120031 ) . Note that 
there is a degeneracy between p* and P, so with a single transit p* must be determined 
by other means in order to derive an estimate of P. p* can be estimated via spectroscopy 
combined with theoretical isochrones, or asteroseismology. The total uncertainty in the 
period is the quadrature sum of the contribution from the estimate of p* and the contribution 
from the K epler photometry (i.e., the uncertainty at fixed p*), 

=f-) 2 + (?) 2 • < 8 > 

P / tot \ P* J V P / Kep 

The contribution from Kepler is dominated by the uncertainty in the ingress/egress time r, 



where we have assumed r<T, and that the number of points taken out of transit is much 
larger than the number taken during the transit (and thus the out-of-transit flux is known 
essentially perfectly). Here Q is approximately equal to the total signal-to-noise ratio of the 
transit, 

Q = (T ph Ty/ 2 6, (10) 
where T p h is the photon collection rate. For Kepler we assume, 

T ph = 7.8 x 10 8 hr- 1 10- a4 ^- 12 \ (11) 



3 As we will show, including limb-darkening and the uncertainty in the out-of-transit flux does not change 
the following discussion qualitatively. 
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and thus Kepler will detect a single transit with a signal-to-noise ratio of 

Q ~ 1300 ( ^V 3/2 ( ^lV 1/6 ( Jjl\ 2 ( _L_) 1/6 i -o.2(v-i2) (12) 

where r p is the radius of the planet. For a Neptune-size planet with these fiducial parameters, 
Q ~ 150. 

Kepler's contribution to the uncertainty in the period is, 

(^:« p — ™— (%)X0\^^py\ (1 a) 

where we have assumed the impact parameter 6 = 0. Due to the strong scaling with r p 
and relatively weak scaling with P, there are essentially two regimes. For r p > i?Nc P , the 
uncertainty in P is dominated by uncertainties in the estimate of p*, which is expected to 
be 0(10%), whereas for r p < R^ep, it is very difficult to estimate the period due to the 
uncertainty in r from the Kepler photometry. 

Once the period is known from the photometry, the mass of the planet, m p , can be 
measured with radial velocity observations. For a circular orbit, and for observations spread 
out over a time that is short compared to P, the radial velocity of the star can be expanded 
about the time of transit, and so approximated by the velocity at the time of transit vq, plus 
a constant acceleration A+, 

v+pavo-A+tt-tc), (14) 

where A+ = 2irK+/P, and 

(2tyG \ 
pM2j m P Sinl ( 15 ) 

is the stellar radial velocity semi-amplitude. Thus, 



G J \2ir J 16CT* \^/6 
where g* is the surface gravity of the star, and we have assumed sin i — 1. 

The uncertainty in the planet mass has contributions from three distinct sources: the 
uncertainty in g+, which can be estimated from spectroscopy, the uncertainty in A+, which 
is derived from RV observations after the transit, and the uncertainties in T, r, and 5, 
which are derived from Kepler photometry. In fact, the uncertainty in r dominates over the 
uncertainties in T and 5. Therefore, we may write, 

fe) 2 -fr) 2+ fe) 2 -(T) 2 <-> 
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For reasonable assumptions and long-period planets (P > 1 yr), we find that the uncertainty 
in A* dominates over the uncertainty in r. 

Assuming that N equally-spaced RV measurements with precision ctrv are taken over a 
time period T tot after the transit, the uncertainty in A* is, 



RV 



A* J AlTl t N' 



- °- 85 (to^t) y-r^) UJ vw) \mJ feJ ,(18) 

where we have assumed N 1. Thus, radial velocity observations combined with Kepler 
photometry can confirm the planetary nature of a Jupiter-sized planet in a relatively short 
time span. An accurate (< 10%) measurement of the mass for a Jupiter-mass planet, or 
even a rough characterization of the mass for a Neptune-mass planet, will require either more 
measurements, or measurements with substantially higher RV precision. Of course, addi- 
tional radial velocity observations over a time span comparable to P will further constrain 
the mass and period. 

In §4] and §0 we derive the above expressions for the uncertainty in the mass and the 
period of the planet using a Fisher information analysis. 



4. Estimating the Uncertainty in P 

4.1. Uncertainties in the Light Curve Observables 

In the absence of significant limb- darkening, a transit light curve can be approximated 
by a trapezoid that is described by the five parameters t c , T, r, 5, and Fq, where Fo is the 
out-of-transit flux. Figure [2] shows this simple trapezoidal model and labels the relevant 
parameters. Mathematically, the flux F as a function of time t is given by, 



F(t) 



Fo, 

Fo - S{t 
F -6, 
F -S(l 
Fo, 



[t e -T/2-T/2])/r, 
[(t-[t c -T/2-r/2])/r]), 



t < ^ 

h<t< t 2 , 
t 2 <t< t 3 , 
t 3 <t< u, 
t 4 < t, 



(19) 



where t\ - are the points of contact. We also define D to be the total duration of the 
observations. This model is fully differentiable and can be used with the Fisher matrix 
formalism to derive exact expressions for the uncertainties in the parameters t c , T, r, 5, 
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and Fq. We note that the definition for 5 we use here differs slightly from the definition we 
adopted in ^J3j Here 6 is the depth of the transit (e.g., in units of flux), rather than the 
fractional depth. These two definitions differ by a factor of F Q such that 8f rac = Sfi ux /F . 
However, when D ^> T (as will be the case for Kepler), the uncertainty in F Q is negligible, 
and so one may define Fo — 1, thus making these two parametrizations equivalent. 

Figure [3] compares th is simple model with the exact model computed using the formalism 
of iMandel fc Agoll (120021 ). in this case for a transit of the Sun by Neptune with a 3. 5- year 
period and an impact parameter of b = 0.2. We see that the trapezoidal model provides an 
excellent approximation to the exact light curve, and as we will show the analytic parameter 
uncertainties are quite accurate. We also show the case of significant limb- darkening as 
expected for a Sun-like star observed in the i?-band (similar to the Kepler bandpass). In this 
case, the match is considerably poorer, but nevertheless we will show the analytic parameter 
unc ertainties we derive assuming the simple trapezoidal model still provide useful estimates 



sec 



Carter et al.l 120081 for a more thorough discussion) . 



We derive exact expressions for the uncertainties in t c , T, r, 5, and F by applying the 
Fisher matrix for malism to the simple light curve model (a detailed explanation of the Fisher 
matrix is given by iGouldl 120031 ) . The full details of the derivation are given in the appendix 
(§EJ). In the case of the long-period transiting planets that will be observed by Kepler, we 
can simplify the full expression by making a couple of assumptions. We assume that the flux 
out of transit, F , is known to infinite precision and that r <C T <C D. The uncertainties in 
the remaining parameters are, 



(7 'J" 



s 




(20) 

(21) 

(22) 
(23) 



We tested the accuracy of these expressions using Monte Carlo simulations. We generate 
1000 light curves for a given set of orbital parameters. We assume that the planet orbits a 
solar-type star in a circular orbit with an impact parameter of 0.2. We repeat the analysis 
for four of the solar-system planets: Jupiter, Saturn, Neptune, and Earth. We assume the 
expected photometric precision for Kepler for a stellar apparent magnitude of V = 12 (a 
representative magnitude for Kepler's stellar sample), a sampling rate of one per 30 min, and 
a total duration for the observations of D — 200 hours. Then we fit for the parameters t c , T, 
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r, 5, and Fq using a down-hill simplex method. The uncertainty in each parameter is taken 
to be the standard deviation in the distribution of the fits for that parameter. A comparison 
of the analytic expressions and the Monte Carlo simulations is shown in Fig. HJ where we 
have plotted the uncertainties as a function of period for a planet with radius 1-RNep crossing 
a star with radius IR® at an impact parameter b = 0.2. We find that the Fisher matrix 
approximation breaks down as the number of points during ingress/egress becomes small. 

We also considered expected un certainties for th e exact uniform source and limb-darkened 



models for the transit light curve. iMandel fc Agoll (120021 ) give the full solution for a tran 



sit involving two spherical bodies and provide code for calculating the transit light curve. 
They also include a limb-darkened solution. These models may be parametrized by the 
same observables described in our simplified model. The full and limb-darkened light curves 
are plotted in Fig. [3] as the dot t ed an d dashed curves, respectively. We use the limb- 



darkening coefficients from I Clare tl (120001 ) closest to observations of the Sun (T e fj- = 5750 K, 
log(g*) = 4.5 cms -1 , [M/H] = 0.0) in the R-band. We apply the Fisher matrix method to 
numerical derivatives of these model light curves to compare the uncertainties for the full 
and limb-darkened solutions with the analytic uncertainties for the trapezoidal model. A 
comparison of the uncertainties from the different models is shown in Fig. |H We find that 
the simplified model is a good approximation for the exact, uniform-source transit model. 
While the comparison is less favorable with the limb- darkening model, the uncertainties do 
not differ by more than a factor of a few. 

We also compared how the uncertainties in the transit observables varied with impact 
parameter b for the three models, because varying b will affect the relative sizes of T and r. 
These comparisons are shown in Fig. [5] for Neptune in a circular orbit around the Sun with 
a 3.5-year period. As can be seen from the figure, variations in b have little impact in the 
uncertainties in the observables for b < 0.8. 



4.2. Uncertainties in Quantities Derived From the Light Curve 



Assuming the planet is in a circul ar orbit, and assuming th e stella r density, p*, is es- 
timated from independent information, ISeager &: Mallen-Ornelad (120031 ) demonstrated that 
the planet period, P, and impact parameter, b, can derived from a single observed transit, 
with sufficiently precise photometry. In fact, provided independent estimates of M* and R± 
are available, then it is also possible to derive the planet radius, r p , semimajor axis, a, and 
instantaneous velocity at the time of transit, t>t r , P - 



The uncertainties in these physical parameters can be calculated through standard error 
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Fig. 2.— The simplified model of a transit light curve showing the points of contact 
(^1,^2,^3,^4), the transit depth (5), the ingress/egress time (r), the FWHM duration of the 
transit (T), the total duration of observations (D), and the time of the transit center (t c ). 
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Fig. 3. — Three different models of a transit light curve. The solid line shows the simplified 
trapezoidal model, the dotted line (barely visible) shows the exact, uniform- source light 
curve, and the dashed line includes the effects of limb -darkening. The e xact uniform-source 
and limb-darkened light curves are calculated using iMandel fc Agoll (120021 ). These light 
curves were generated for a Neptune analog orbiting a solar analog at a period of P = 3.5 yrs 
and an impact parameter b = 0.2. Typical error bars for Kepler assuming V = 12.0 and 
30-min sampling (solid) or 20-min sampling (dotted) are indicated. 
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Fig. 4. — Uncertainty in observables as a function of period. The uncertainties in t c , T, 
t, and 5 are plotted versus period for Neptune orbiting the Sun with an impact parameter 
of 0.2. The different line styles indicate the trapezoidal transit model (solid line), the full 
solution (dotted line), and a model including the effects of limb-darkening (dashed line). 
Monte-Carlo simulations of the trapezoidal model using 30-minute sampling are shown as 
crosses. The gray crosses show simulations for 20-minute sampling. The dotted line in the 
bottom panel is not visible because it is nearly identical to the solid line. The vertical dotted 
lines mark where the ingress/egress duration r is equal to the sampling rate dt, as well as 
where r = 2dt for 30-minute sampling. 
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Fig. 5. — Uncertainty in observables as a function of impact parameter. The uncertainties 
in t c , T, t, and S are plotted versus impact parameter for Neptune orbiting the Sun with a 
period of 3.5 years. The line styles are the same as for Fig. HI As in Fig. m the dotted line 
in the bottom panel is not visible because it is nearly identical to the solid line. 
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propagation. Equations for these quantities and their uncertainties are given in Eqs. [23]- 
[281 (Note that we calculate b 2 instead of b because numerical methods have difficulty when 
Ts/5/r ~ 1). In each case below, the second approximate uncertainty equation is calculated 
by substituting the variances and covariances of the observables r, T, and S, and simplifying 
under the assumption that r <C T <C D. 
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Of course, the quantities i?*, gv, and p* and their uncertainties must be estimated from 
some other external source of information, such as spectroscopy or theoretical isochrones. 
We can generally expect that the fractional uncertainties on these quantities will be of 
order 10%. For long-period (P > L = 3.5 yrs), Jupiter-sized planets, Q > 1, and thus 
the variances and covariances of the transit observables will be small in comparison to the 
expected uncertainties on -R*,p*, and g+. 

The fractional uncertainties in these derived parameters are plotted in Fig. [61 We 
compared the analytic uncertainties given in Eqs. [2414281 to the Monte Carlo simulations 
described in Sec. 14.11 The fractional uncertainty in each parameter from the Monte Carlo 
simulations is one-half the range of the middle 68% of the data divided by the median. Notice 
that as r decreases from 2dt to dt, the simulations become increasingly disparate from the 
theoretical expectations. ov p /r p does not show this behavior because it does not depend on 
r. In this regime, as r decreases, it becomes increasingly probable that only one point will 
be taken during ingress leaving r relatively unconstrained and increasing its uncertainty and 
the uncertainty of quantities that depend on it. 

We now consider in detail the uncertainty in the estimated period. Given that Q oc 
T 1//2 <5, we have that Q~ 2 (T/r) oc (5 2 r)~ 1 . Since r oc r p P l l z and 5 oc r 2 ,, the contribution to 
the uncertainty in the period due to the Kepler photometry is a much stronger function of 
the planet's radius than its period, Q~ 2 (27T/2r) oc r~ 5 p -1 / 3 . As a result, for a fixed stellar 
radius i?*, the ability to accurately estimate P depends almost entirely on r p . Furthermore, 
because the scaling with r p is so strong, there are essentially two distinct regimes: for large r p , 
the uncertainty in P is dominated by the uncertainty in p*, whereas for small r p , the period 
cannot be estimated. The boundary between these two regimes is where the contribution to 
the uncertainty in the period from p* is equal to the contribution from the K epler photometry, 
i.e. where the two terms in the brackets in Eq. [28] are equal. This occurs at a radius of, 



p,crit — ^Nep 



10 0.4(V-12) 



(Tjp*Y 2 ( R*V ( M*\ 1/3 /3.5yrs^ 1/3 l 



0.1 ) \R Q J \M q J \ P 



(29) 



assuming 6 = 0. Note that r P)Cr i t is a weak or extremely weak function of all of the parameters 
except R*. 

Fig. [7| illustrates this point by showing contours of constant uncertainty in the period 
as a function of the radius and period of the planet. The calculations are for systems with 
R± = Rq, b = 0.2, and V = 12.0. This figure shows that for planets larger than the radius of 
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Neptune, the uncertainty in the period will be dominated by the uncertainty in the inferred 
stellar density, whereas for planets much smaller than Neptune, an accurate estimate of the 
period from the Kepler light curve will be impossible. 

As shown above, the uncertainty in P is dominated by the uncertainty in p± and r. Our 
ability to determine r, and its uncertainty, depends both on its length and how many points 
we have during the ingress or egress. The length of r depends on f tr p (which is a proxy for 
P in the circular case), b, and r p . These properties are intrinsic to the system, but may be 
derived from the observables. We can explore how the contours shown in Fig. [7J vary with 
impact parameter. Figure [5] shows that a system with a larger impact parameter will have 
a smaller fractional uncertainty in r. Thus, a system with a larger impact parameter would 
have a smaller fractional uncertainty in P at fixed period and planet radius. 

The other effect that influences the contours is the sampling of r. As the number of 
points taken during the ingress/egress becomes small, the fractional uncertainty in r, and 
hence in P, increases. In particular, if only one point is taken during the ingress, then the 
duration is relatively unconstrained. The probability of taking only one point during the 
ingress increases linearly from to 1 as r decreases from 2dt to dt. In Fig. El the Monte 
Carlo simulations for the fractional uncertainty in P diverge from the simple model over this 
range. We indicate this region by the shaded portion of Fig. [7J Where r = dt, the fractional 
uncertainty in P can be several times that predicted by the theoretical calculations, but 
these uncertainties converge as the shading gets lighter towards r = 2dt. Below the line 
where r = dt, the uncertainty in r increases rapidly, but fortunately, in much of this regime, 
Kepler will observe multiple transits, and this analysis will be unnecessary. As shown in 
the figure, a faster sampling rate, such as 1 per 20 min, significantly expands the parameter 
space over which our theoretical uncertainties are valid. 



The mass of the planet comes from sampling the stellar radial velocity curve soon after 
the transit is observed. Near the time of transit we can expand the stellar radial velocity, 



where Vq is the systemic velocity, is the stellar radial velocity semi-amplitude, and = 
2ttK*/P is the stellar acceleration. Because the planet is known to transit and has a long 
period, we assume sinz = 1. A Fisher matrix analysis of the linear form of Eq. [30] gives the 



5. Estimating the Uncertainty in m p 
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Fig. 6. — Uncertainty in derived quantities as a function of period. The uncertainties in r p , 
a, P, i>tr,p! an d b 2 are plotted versus period for Neptune orbiting the Sun with an impact 
parameter of 0.2. The line styles are the same as for Fig. HI The vertical dotted lines mark 
where the ingress/egress duration r is equal to the sampling rate dt, as well as where r = 2dt 
for 30-minute sampling. 
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Fig. 7. — Contours of constant fractional uncertainty in the period as a function of the 
period and the planet radius. The model is for a solar-type star with a V magnitude of 12.0 
and an impact parameter for the system of 0.2. The dashed line indicates how the fractional 
uncertainty in P changes with V magnitude; it shows the 0.20 contour for V = 14.0 and 
b = 0.2. The diagonal black dotted line represents the boundary below which the assumptions 
of Fisher matrix break down (the ingress/egress time r is roughly equal to the sampling rate 
of 1/30 min). Thus, the contours below this boundary are shown as dash-dotted lines. The 
region between r = 2dt and r = dt is shaded to indicate the increasingly probability of 
obtaining only one point during ingress or egress leading to an uncertainty in P that is 
larger than theoretical expectations. The diagonal gray dotted line indicates the r — dt 
boundary for 20-min sampling. The vertical dotted line shows the mission lifetime of Kepler 
(L = 3.5 yrs). Solar system planets are indicated. 
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estimated uncertainty in A* to be 



2 _ I^rv ,,>,, 

1 tot JV 

where ctrv is the radial velocity precision, T tot is the total time span of the radial velocity 
observations, and N is the number of observations, and we have assumed N ^> 1 and that 
the observations are evenly spaced in time. The details of this derivation are given in the 
appendix (£JB]). 

Equations for the mass of the planet and the uncertainty are 
1 , . fTrV 



™v = g:A (32) 



2 2 



2 "^2 ^2 ^2 2 ^2 ^2 ^2 



4 2 1 2 ^2 /^24T 



The approximate expression for the uncertainty in the mass is taken in the limit r <C T <C D. 
Furthermore, the Q~ 2 (24T/r) term can be neglected in most cases because, as we showed 
for the corresponding term in uncertainty in the period, it will be very small for planets with 
radii larger than that of Neptune. The scaling for the (cr^/A*) 2 term is given in Eq. [TH 

Contours of constant fractional uncertainty in the mass of the planet are shown in Fig. 
[8] as a function of P and m p , assuming N = 20 radial velocity measurements are taken 
over T to t = 3 mos with a precision ctrv = 10 ms -1 . It shows that the planet mass can be 
estimated to within a factor of two over this time period, thus establishing the planetary 
nature of the transiting object. By doubling the length of observations, one can put stronger 
constraints on the mass of the planet. 

There are two points to bear in mind when applying this estimate for the uncertainty in 
the mass of the planet. First, for r p < -RNeptune, the Q~ 2 (24T/t) term is no longer small. The 
second consideration is that as time progresses away from the time of transit, the straight 
line approximation to the radial velocity curve will break down. In that case, the period will 
begin to be constrained by the radial velocity curve itself, and the uncertainties should be 
calculated from a Fisher matrix analysis of the full expression for the radial velocity curve 
with three parameters: Vq, K+, and P. 
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Fig. 8. — Contours of constant uncertainty in m p as a function of m p and P. The solid lines 
show the result for 20 radial velocity measurements with precision of 10 m s _1 taken over a 
period of 3 months after the transit. The dashed line shows the contour for <J mp /m p = 0.50 
for 40 observations taken over 6 months. The dotted line indicates the mission lifetime of 
3.5 years. The positions of Jupiter and Saturn are indicated. 
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6. Eccentricity 

The results presented in the previous sections assumed circular orbits. Given that 
the average eccentricity e of planets with P > 1 yr is e ~ 0.3, this is not necessarily a 
good assumption. In this section, we assess the effects of non-zero eccentricities on the 
ability to characterize single-transit event s detected by Kep ler. We note that our discussion 



has some commonality with the study of IFord et al.l (120081 ). who discuss the possibility of 



characterizing the orbital eccentricities of transiting planets with photometric observations. 
However, our study addresses this topic from a very different perspective. 

In §3] we demonstrated that, under the assumption that e = 0, the planet period P and 
mass m p can be inferred from two observables: the velocity of the planet at the time of transit 
v tr, and the projected acceleration of the star A*. In the case of a non-zero eccentricity, there 
are two additional unknown parameters: the eccentricity e and the argument of periastron of 
the planet u p . Thus, with only two observables (v tT and A+) and four unknowns (P, m p , e, u p ), 
it is not possible to obtain a unique solution. 

Although a unique solution to the planet parameters does not exist using the observable 
parameters alone, it may nevertheless be possible to obtain an interesting constraint on P 
and/or m p by adopting reasonable priors on e and u p . From the transit observables, one can 
estimate the velocity of the planet at the time of transit (Eq. [5]). In the case of e ^ 0, this 
is given by, 

_ 2na 1 + e sin u p 

ftr.p - p q _ e 2y/ 2 ' \ 66 > 



where here and throughout this section we assume sinz = 1. Solving for P, 



2nGM± 

V tr,p 



1 + e sin u p 

- P.2U/2 



(34) 



For a fixed value of e, the inferred period relative to the assumption of a circular orbit has 



extremes due to the unknown value of u p of 



P \ 1 

/ mini max \ 

Taking a typical value for the eccentricity of e = 0.3, this gives a range of inferred periods of 
0.4 — 2.5 relative to the assumption of a circular orbit. In fact, because the probabili t y that 



a pla net with a given a transits its parent star depends on e and u p (jBarned 120071 ; iBurke 



20081 ). a proper Bayesian estimate, accounting for these selection effects, might reduce the 



range of inferred values for P significantly. 
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What can be learned from RV observations immediately after the transit? Expanding 
the stellar projected velocity around the time of transit, we can write, 

v+ = v + eK* cosu p - A*(t - t c ) + -J*(t - t c f. (36) 
The projected acceleration of the star at the time of transit is, 

A= ^(i±£^y. ,37) 

cr \ 1 — e z J 
This can be combined with Vt r , p to derive the mass of the planet, 

GM 2 

m p = — j-J^A^l + esinu p ) 2 (38) 

V tr,p 

Thus, for a fixed eccentricity, the inferred mass relative to the assumption of a circular orbit 
has extremes of, 

— 2 =(l±e) 2 , (39) 

777 / 

P / min/max 

which for e = 0.3 yields a range of 0.5 — 1.7, which is generally smaller than the contribution 
to the uncertainty in m p due to the measurement uncertainty in A+ for our fiducial case of 
m p = M Jup and P = 3.5 yrs (see Eq. [18]) . 

We may also consider what can be learned if it is possible measure the curvature of the 
stellar radial velocity variations immediately after transit. The projected stellar jerk is given 
by, 

47rGm p (1 + e sinci;p) 3 e coscUp , . 

J * = ~Ta~ 2 (1 - e 2 ) 7 /2 • (4U) 

This can be combined with t> tr and A* to provide an independent constraint on a combination 
of the eccentricity and argument of periastron, 

ecosa;p _ GM+J* 
(1 + esinu p ) 2 ~ 2A+vl ' ^ ' 

Unfortunately, it will be quite difficult to measure J*. Using the same Fisher formalism as 
we used to estimate the uncertainty in (§HD, assuming N evenly-spaced RV measurements 
with precision ctrv, taken over time span T tot after the transit, the uncertainty in J* is, 

2 ^ 720^ 



< - (42) 



NT 4 
1 v 1 tot 
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where we have assumed N 3> 1. For our fiducial parameters, the fractional uncertainty in 
J* is, 

fz^Y ~ 710 ( gRv v p mos Y w f^Y f^Y /3 f p Y° /s ( ecosuj p y 2 

\J*J UOms- 1 / V T tot j V ™ p / Wo/ \3-5yrsJ V 0.3 / 

(43) 

where we have approximated ecosu;p(i + e sinc<jp) 3 (l — e 2 )~ 7 / 2 ~ ecoscj p . We conclude that, 
in order to obtain interesting constraints on the eccentricity of long-period planets detected 
by Kepler, higher-precision RV measurements taken over a baseline comparable to ~ P will 
be needed. 



7. Summary 

The discovery of long-period transiting planets along with subsequent follow-up obser- 
vations would greatly enhance our understanding of the physics of planetary atmospheres 
and interiors. Such planets would allow us to gather constraints in a regime of parameter 
space currently only occupied by our own giant planets, namely planets whose energy bud- 
gets are dominated by their residual internal heat, rather than by stellar insolation. These 
constraints, in turn, might provide new insights into planet formation. In this paper, we 
demonstrated that it will be possible to detect and characterize such long-period planets 
using observations of single transits by the Kepler satellite, combined with precise radial 
velocity measurements taken immediately after the transit. Indeed, these results can be gen- 
eralized to any transiting planet survey using the scaling relations we provide, and it may 
be particularly int eresting to ap ply them to the COnvection, Rotation & planetary Transits 



(CoRoT) mission flBagrinll200 



app 
)J). 



We calculated that Kepler will see a few long-period, single transit events and showed 
that, for circular orbits, the period of the system can be derived from the Kepler light curve 
of a single event. We derived an expression for the uncertainty in this period and showed that 
it is dominated by the uncertainty in the stellar density derived from spectroscopy (which we 
assume is ~ 10%) for planets with radii larger than the radius of Neptune, rather than being 
dependent on the properties of the transit itself. This method can also be applied to planets 
that Kepler will observe more than once, so that the second time a planet is expected to 
transit Kepler can make a selective improvement to its time sampling to better characterize 
the transit. 

We have also shown how the mass of the planet can be constrained by acquiring precise 
(~ 10 m s _1 ) radial velocity measurements beginning shortly after the transit occurs. We 
have shown that 20 measurements over 3 months can measure the mass of a Jupiter-sized 
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object to within a factor of a few and that extending those observations to 40 measurements 
over 6 months significantly reduces the uncertainty in the mass. Knowing the mass to 
within a factor of a few in such a short time can distinguish between brown dwarfs and 
planets rapidly and allow us to maximize our use of radial velocity resources. 

We explored the effect of eccentricity on the ability to estimate the planet mass and 
period. Allowing for a non-zero eccentricity adds two additional parameters, and as a result 
it is not possible to obtain a unique solution for the planet mass, period, eccentricity, and 
argument of periastron. However, by adopting a reasonable prior on the eccentricity of the 
planet, the period and mass of the planet can still be estimated to within a factor of a few. 
Detailed characterization of the planet properties will require precision RV measurements 
obtained over a duration comparable to the period of the planet. 

Thus, in the interest of "getting the most for your money," we have shown that the 
sensitivity of Kepler extends to planets with periods beyond its nominal mission lifetime. 
With the launch of the Kepler satellite, we are poised to discover and characterize several 
long-period transiting systems, provided that we are prepared to look for them. 

We are grateful to Josh Carter, Jason Eastman, Eric Ford, Andy Gould, Matt Holman, 
Yoram Lithwick, Josh Winn, and Andrew Youdin for helpful discussions. We thank the 
referee Ron Gilliland for constructive comments. JCY is supported by a Dean's Graduate 
Enrichment Fellowship from The Ohio State University. 
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8. Appendix: Detailed Derivations 

8.1. Derivation of the Uncertainties in the Light Curve Observables 

The Fisher matrix formalism is a simple way to estimate the uncertainties in the param- 
eters, a, of a model F(x), that is being fit to a series of measurements x k with measurement 
uncertainties a k . The covariance of with ctj is given by the element Qj of the covariance 
matrix c, where c = b" 1 and the entries of b are given by, 



^ dF{x k ) dF(x k ) 1 

k=l 



where is the number of data points. In the limit of infinite sampling (N^ 
precision, a k = cr, 

" D dF(x) dF(x) 



b i:j 



(44) 

oo) and fixed 
(45) 



Da 2 J don dctj 

where the interval of interest is given by x = [0, D] , and in this case D is the total duration 
of observations. Thus, if the partial derivatives of a model with respect to its param eters 
are known, then the uncertainties in those parameters can be estimated (jGouldll2003l ). 



For the simplified transit model described in §4.11 the observable parameters are a 
[t c , T, r, S, Fq]. For a sampling rate 7, the b matrix is, 
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The covariance matrix is, 
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(47) 

where tf = T — r is the duration of the flat part of the eclipse and t out = D — t — T is the 
time spent out of eclipse. Thus, the uncertainty on the ingress/egress time, a T , is given by 



0V 



a -t(5Dt - At 2 - 6DT + 6T 2 ) 



(48) 



Define Q = y/ r yT{8/a). Then, in the limit r<T, the uncertainties in the observable 
parameters are given by, 



1± 

5 




= , since in general, 5 <C F . 



(49) 

(50) 

(51) 

(52) 
(53) 



Note that Q is approximately the total signal-to- noise ratio of the transit. Assuming that 
the photometric uncertainties are limited by photon noise, we have that 7/a 2 = r p h, where 
r p h is the photon collection rate. This recovers the expression for Q in §3j 

A more detailed ana l ysis o f the variances and covariances of the transit observables can 
be found in ICarter et al.l ( 2008 ). 



8.2. Derivation of the Uncertainties from the Radial Velocity Curve 



A similar analysis can be done for the radial velocity curve, except that we use Eq. HH 
and consider discrete observations. The radial velocity for a circular orbit is 
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Expanding about t c gives 



v ir = v - sin 



2tt 



2tt 
~P 



(t - t c ) 



v* = v - K*—(t - t c ) =v - A+(t - t c ), 

where A+ is the stellar projected acceleration at the time of the transit, 

2irK+ 



A, 



(54) 
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(56) 



Consider N measurements of u*, each with precision ctrv, taken at times Fitting 
these measurements to the linear model in Eq. [551 we can estimate the uncertainties in the 
parameters vq and A* using the Fisher matrix formalism, 
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The covariance matrix is the inverse of b, 
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The uncertainty in the projected stellar acceleration is, 



'RV 



N(N — 1)(JV + l)(At) 2 ' 



(62) 



In the limit as iV — > oo, the covariance matrix reduces to, 



(T 



C = 



RV 



At 2 



4At 2 6At 

N N 2 

6At 12_ 

N 2 N 3 



(63) 



Defining the total length of observations as T tot = AT At, we can also write, 



a 



RV 



T 2 

1 tot 



A Ttot a Ttot 

N u N 
a Ttot 12 
u N N 



(64) 



and so the uncertainty A* in the limit of N — > oo becomes, 



12o|v 

T t 2 ot iv • 



(65) 



